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Within the scope of a LRS Bianchi type-I cosmological model we study the role of the 
nonlinear spinor field in the evolution of the Universe. In doing so we consider a polynomial 
type of nonlinearity that describes different stages of the evolution. Finally we also use the 
observational data to fix the problem parameters that match best with the real picture of 
the evolution. The assessment of the age of the Universe in case of the soft beginning of 
expansion (initial speed of expansion in a point of singularity is equal to zero) the age was 
found 15 billion years, whereas in case of the hard beginning (nontrivial initial speed) it was 
found that the Universe is 13.7 billion years old. 
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I. INTRODUCTION 


The discovery and further confirmation of the existence of the accelerated mode of expansion 
of the present day Universe [T, ^ lead cosmologists to construct new theories able to explain 
this new cosmological findings, T rough cosmological constant, quintessence, Chaplygin gas etc. 
are the prime candidates, thank? fi( a number of remarkable works [3-171. recently many authors 
considered the spinor fields as a possible alternative to these mode 5. And it is because of the 
spinor fields’ ability to simulate different type of source fields frorrLf ^ Totic matter to phantom 
matter and Chaplygin gas [18-221. Moreover it was found that a nonlinear spinor field can also 
(i) generate singularity-fre J Uiherse [^7, 9, 1^; (ii) accelerate the isotropization process of 
the initially anisotropic spa|ce|iprq [7, 9, L2| ai(l(Tii) give rise to a late time accelerated mode of 
expansion [11. 13-151. | 

Some re cert sn dijs show that tpy |n|o|n2|[iagonal components of the energy-momentum tensor 


of the spingr field 


qanj play significant ro 


_ e on the geometry of spacetime as well as on the compo¬ 
nents of the spinor field itself, namely on the spinor field nonlinearity [23--301 . In those papers it 


was shown that depending on the specificity of metric in some cases tl e f rtin 
and the mass terms vanish all together, whereas in other cases both th(; rqi^ss 


Dr field nonlinearity 
term and nonlinear 

term do not only disappear but also play a crucial role in the evolution of the Universe. In a re¬ 
cent paper [301 it was found that within the scope of a LRS Bianchi type-I spacetime the spinor 
field nonlin DMty depending on the sign of self-coupling constant allows either an expansion with 
acceleratior _oj an oscillatory mode of evolution of the Universe. In this note we develop that the¬ 
oretical work further and fix the values of the problem parameters using the recent observational 
data. 


II. BASIC EQUATIONS 


The LRS Bianchi type-I (BI) model is the ordinary Bianchi type-I model with two of the three 
metric functions being equal to each other and can be given by 


ds^ = dt^ — af [dx^ + dy^] — dz^, 


( 2 . 1 ) 


with ai and aj, being the functions of time only. 

Keeping this in mind the symmetry between y/ and y/ we choose the symmetrized Lagrangian 
[311 for the spinor field as : 


_ i 

L=- 




mspy/y/-F, 


( 2 . 2 ) 


where the nonlinear term F describes the self-interaction of a spinor field and can be presented as 
some arbitrary functions of invariants generated from the real bilinear forms of a spinor field. We 
consider F = F{K), with K = {I, J, I+ J, I — J}. It can be shown that such a choice describes the 
nonlinearity in its most general form. 

Varying (12.21) with respect to yr{y^) one finds the spinor field equations: 

lyf^V^y/ - mspy/ - 2 Fk{SKi + iPKj)^)y/ = 0, (2.3a) 

ly^yrf +m,pyf + 2FKnSKj + iPKjy^)=0. (2.3b) 

Here we denote F^ = dF/dK, Kj = dK/dl and Kj = dK/dJ. 

The energy-momentum tensor of the spinor field is given by 


'T'r _ 

J- II — 


= -gP^ 
4« 


¥Yvyfi¥-yfi¥rv¥-^v¥rfi¥] - 


-sp 


(2.4) 
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where L^p in view of (12.31) can be rewritten as 

= -m,p\frxir-F{K) 

= 2Fk{IKi + JKj) -F = 2 KFk -F{K). (2.5) 

It can be shown that the spinor field in this case possesses nontrivial non-diagonal components 
of the energy-momentum tensor. On account of that the system of Einstein equations can be 
written as 1301 


+ — + —— = k(F(K) - 2KFk), 

(2.6a) 

ai as ai 

2^ + ^ = k(F(K)-2KFk), 

(2.6b) 

ai aj 

4 + 2—— — K(m^„S + F(K)), 

(2.6c) 

ci 2 ci\ 

V«3 a\J 

(2.6d) 

V“1 “3/ 

(2.6e) 


where y/ are the components of the pseudovector. 


III. SOLUTION TO THE FIELD EQUATIONS 

In this section we solve the eq uation s obt ained in the previos section. The off-diagonal com¬ 
ponents of the Einstein equations (I2.6dl) and (I2.6el) impose the following restrictions either on the 
components of the spinor field or on the metric functions: 

A^ = 0, a 1 = 0 , 

^ ^ - 0 
as ai 

Erom (I3.1bl) we dully find as = qoai with qo being some constant. In this case the system can 
be described by a ERW model from the very beginning. Hence we don’t consider this case here. 
We will do it in so me of o ur forthcoming papers on ERW model. 

So we cons ider (I3.1al) when the restriction is imposed on the components of the spinor field. 
Subtraction of (I2.6bh from (I2.6al) gives 


(3.1a) 

(3.1b) 


+ ) - 0 
as ai ai \as ai ’ 


(3.2) 


that leads to [71 


dt 


4(1 =Diy^/^exp (Xi / y], as = (l/L>?)y^/^exp f-2Ai j y 


(3.3) 
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with Di and Xi being the integration constants. Thus we see that the metric functions can be 
expressed in terms of V. 

Ou r next step will be to define V. Combining the diagonal Eistein equations (I2.6al) . (I2.6bl) and 
(I2.6cl) in a certain way for V we find [71 


V 



S + 2{F{K)-KFk))V. 


(3.4) 


Now order to solve (13.41) we have to know the relation between the spinor and the gravitational 
fields. Using the equations 


5o + ^a|] = 0, (3.5a) 

/’o-4>A[j = 0, (3.5b) 

where we denote Sq = SV, Pq = PV it can be show that 

y 2 

K = {IJJ + J,I-J}. (3.6) 

The relation (13.61) holds for ^ = {J, I+J, I — J} only for massless spinor field, while for .ST = 7 it 
holds both for massless and massive spinor field. In case of K = I + J one can write S = sin (Vq/V) 
and P = cos(Vo/3^)> whereas for K = I — J one can write S = cosh(V()/U) and P = sinh(Vo/l^). 
In what follows, we will consider the case for K = I, setting F = as in this 

case further setting spinor mass m^p = 0 we can revive the results for other cases. 

Then inserting F = = 'Lk^kS^"’' into (13.41) and taking into account that in this case 

S = Vq/V we find 




m 


spyo+2£A,(i-n,)yo"”^y 
k 




with the solution in quadrature 


dV 


3k 


m 


sp 


vbv+E^A^yo^”'y2(i-«^) fc 


— t + to, 


(3.7) 


(3.8) 


with C and to being some arbitrary constants. 

In what follows we solve the equations for V , i.e., (13.71) numerically. But before doing that we 
write the solution to the spinor field equations exp licitly. 

The solutions to the spinor field equations (12.31) in this case can be presented as 

¥l 2 (f) = ’ ¥3,4(0 = ’ (^-^0 

with Cl, C 2 , C 3 , C 4 being the integration constants and related to yo as 

^ Cl + C2C2 — C3C3 — C4C4 — KO* 

Thus we see that the metric functions, the components of spinor field as well as the invariants 
constructed from metric functions and spinor fields are some inverse functions of V of some de¬ 
gree. Hence at any spacetime point where y = 0 it is a singular point. So we con side r the initial 
value of y (0) is small but non-zero. As a result for the nonlinear term to prevail in (13.71) we should 
have 1 — Ink < i.e.> > 1/2, whereas for an expanding Universe when y —)■ 00 as t —)■ 00 one 

should have 1 — 2ni_ > 0, i.e., < 1/2. As is seen from (13.71) . n]^ = 1/2 leads to a term that can 
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be added to the mass term. So without losing the generality we ean eonsider no = 1/2, ni =0 and 
n2 = 2 . 

In this ease we obtain 


3 fc 

y = d>i(y), $i(y) = — [(m^p + Ao) Vq+2XiV -2X2V^v~^]. 


Equation (13.101) allows the first integral 


y = d> 2 (y), 4>2(y) = [(msp+Ao) yoy+Aiy 2 +A2yo4y-2+c]. 

The solution to the equation (13.101) ean be written in quadrature as follows 

dV 


4>2(y) 


— t + to- 


(3.10) 


(3.11) 


(3.12) 


To solve the (13.101) we should ehoose the problem parameters Vq, mgp, k, C, A^ as well as the 
initial value of V (0) in sueh as way that does not leads to 

(msp + Ao) Toy + Aiy2 + X 2 V^V-^ + C<0. 

In a reeent paper [301 we eonsidered the ease with Aq = Ai = A 2 = A. It was shown that in ease 
of positive A we have ai aeeelerated mode of expansion of the Universe, while for negative A we 
have oseillatory solution. 

In this paper we do not perform numerieal analysis to obtain different ty pe of solutions for 
different values of problem parameters. In what follows we study the equation (13.101) numerieally 
to find the problem parameters that fit best with the observational data. 


IV. COMPARISON WITH OBSERVATIONS 


In what follows we numerieally solve the equation (13.71) . Let us rewrite this equation as follows: 

3 

y=A+^5,(l-nfc)yi-2«.. (4.1) 

k=l 

Here A = ^ (msp + Aq) yo, i-e., the eonstant term, whereas = 2XkVQ"''. 

Defining Hubble parameter and red-shift as follows 


(4.2a) 
(4.2b) 

we rewrite (14.11) in the form 

i = -i W - 1)''"*-^ 

(4.3) 

d 1 

[z+IAVh 

with initial values t(0) = 0, H{0) = Hq. 


3y a 
1 \ 1/3 
Z+l=(^ 
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The foregoing equation allows the following solution 


HM = 




^A[z.+ 1)^ + 5 £ Bt(z+ l)'S”i +Ch(z+ 1)«, 


(4.4) 


k=l 


Ch is eonstant of integration. 

The numerieal values of the parameters sueh as spinor mass, power of nonlinearity ete. are 
determined by eomparing the solutions to astrophysieal observations exploiting the maximum 
likelihood method by minimizing the funetional 




(4.5) 


where Zi, Hi and Ohi were taken from the tables given in [32-341. 


From (I4.2bl) it follows that °° leads to a{t) —)■ 0, i 
For an expanding Universe it means z —)■ oo at the time of 
(soft origin). From d = ^ 


there oeeurs a spaee-time singularity. 
^IBang. Let us assume that d\^^oo — 0 
From d = —)■ 0 we see that it may happen it z inereases faster than H (not 


neeessarily H = 0). Then from (14.41) we ean eonelude that A = 0, Ch = 0 and < \ /3. The 
equality A = 0 states that in this ease the spinor mass is either zero or eompensated by some 
eonstant. Analogieal results were found in the papers where dark energy was simulated by spinor 
field [271 . 

The|alue of is searehed using the method of random walk with the seleetion of sueeessful 
steps.L^t the beginning =0.1 and At eaeh m-th step we replaee = 

+0.01(^ using standard normal distribution of random variable ^ G N{Q, 1), but under the 
eondition 

To determine A and at eaeh m-th step for fixed njc the funetional (14.51) is approximately 
substituted by 


what is quadratie to eoeffieients to be determined and the quest of the minimum is redueed to the 
solution of linear equations. In Fig. [T] we have plotted the dynamies of the proeess of finding the 
minimum X- 


. 

5 10 15 20 25 30 35 


FIG. 1. Dynamics of the process of finding the minimum x^ 
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The adjusted values of parameters are 

m =0.1961075852, = 8.033381 ■ 10^ = 

^2 = 0.1922624717, B 2 = -1.073812-10^ = 3^fA2yo^”^ (4-V) 

^3 = 0.1805648222, = 2.709600-10^ = 

At present epoeh model time t = 0, red-shift z = 0 and Hubble parameter Hq = 13.5. 

It is known that z = Ht for small z, so eomparing at red-shift z = 0.001 the astronomieal 
time t“( 0 . 001 ) = O.OOI/T/q with model time t(O.OOl) we ean obtain the age of Universe T = 
[t(oo)/t(0.001)]t"(0.001) = 15.0-10^ years. The eonfidenee level is found to be CL = 0.92. 

In ease of A 7 ^ 0 (hard origin), the adjusted values of parameters are 

A = -39093.17887 = ^ (m^p -f Aq) Vq, 
ni =0.3318171453, 5i = 10.70355153 -10^ = , 

^2 = 0.3122973942, ^2 =-21.81174544-10^ = 3^fA 2 yo^"^ ' 

R 3 = 0.2939425415, B 3 = 11.23402881 -10^ = SkX^Vq'^L 

At present epoeh model time t = 0, red-shift z = 0 and Hubble parameter T/q = 72.8, we ean obtain 
the age of Universe T = 13.7 -10^ years. The eonfidenee level is found to be CL = 0.93. 

In Figs. |2]and|3l we have plotted the H{z) /(1 -t-z) data (32 points) and model predietion (line for 
best-fit model) as a funetion of red-shift (logarithmie seale) for hard and soft origins, respeetively. 
In Figs. |4]and[5]//(z)/(l -t-z) model predietion as a funetion of time t (years) has been drawn 
for hard and soft origins, respeetively. Figs. | 6 ] and |7] we have demonstrated the evolution of 
volume seale V (logarithmie seale) model predietion as a funetion of time t (years) for hard and 
soft origins, respeetively. 



FIG. 2. //(z)/(l -t-z) data (32 points) and model 
prediction (line for best-fit modef) as a function 
of red-shift (logarithmic scale) for hard origin 



FIG. 3. H(z)/(l+z) data (32 points) and model 
prediction (line for best-fit modef) as a function 
of red-shift (logarithmic scale) for soft origin 
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FIG. 4. H{z)/{l+z) model prediction as a func¬ 
tion of time t (years) for hard origin 





FIG. 6. V (logarithmic scale) model prediction 
as a function of time t (years) for hard origin 


80-1 
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FIG. 5. H{z)/{'i-+z) model prediction as a func¬ 
tion of time t (years) for soft origin 



FIG. 7. V (logarithmic scale) model prediction 
as a function of time t (years) for soft origin 


V. CONCLUSION 

Within the scope of LRS Bianchi type-I model we have studied the role of spinor field in the 
evolution of the Universe. Since the expression for metric functions, components of the spinor field 
and invariants constructed from these quantities are in functional dependence on volume scale V 
it is important to find and solve the equation for V. In this report we have solved the equation 
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in question and in doing so we have exploited the astronomical data available to fix the problem 
parameters such as coupling constants etc. that provides the best correspondence with observation. 

The assessment of the age of the Universe in case of the soft beginning of expansion (initial 
speed of expansion in a point of singularity is equal to zero) the age was found 15 billion years, 
whereas in case of the hard beginning (nontrivial initial speed) it was found that the Universe is 
13.7 billion years old. 
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